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The  seasonal  ocean  circulation  and  thermal  structure  in  the  South  China  Sea  (SCS)  were  studied 
numerically  using  the  Princeton  Ocean  Model  (POM)  with  20  km  horizontal  resolution  and  23  sigma 
levels  conforming  to  a  realistic  bottom  topography.  A  sixteen  month  control  run  was  performed  using 
climatological  monthly  mean  wind  stresses,  restoring  type  surface  salt  and  heat,  and  observational 
oceanic  inflow/ outflow  at  the  open  boundaries.  The  seasonally  averaged  effects  of  isolated  forcing 
terms  are  presented  and  analyzed  from  the  following  experiments:  1)  non-linear  dynamic  effects  re¬ 
moved,  2)  wind  effects  removed,  3)  open  boundary  inflow/  outflow  set  to  zero.  This  procedure  allowed 
analysis  of  the  contribution  of  individual  parameters  to  the  general  hydrology  and  specific  features  of 
the  SCS,  for  example,  coastal  jets,  mesoscale  topographic  gyres,  and  counter  currents.  Our  results  show 
that  the  POM  model  has  capability  of  simulating  seasonal  variations  of  the  SCS  circulation  and 
thermohaline  structure.  The  simulated  SCS  surface  circulation  is  generally  anticyclonic  (cyclonic)  during 
summer  (winter)  monsoon  period  with  a  strong  western  boundary  current,  with  a  mean  maximum  speed 
of  0.5  m/  s  (0.95  m/  s),  mean  volume  transport  of  5.5  Sv  (10.6  Sv),  extending  to  a  depth  of  around  200 
m  (500  m).  During  summer,  the  western  boundary  current  splits  and  partially  leaves  the  coast;  the 
bifurcation  point  is  at  I4°N  in  May,  and  shifts  south  to  10°N  in  July.  Besides,  a  mesoscale  eddy  in  the 
Sunda  Shelf  (Natuna  Island  Eddy)  was  also  simulated.  This  eddy  is  anticyclonic  (cyclonic)  during  sum¬ 
mer  (winter)  monsoon  with  maximum  swirl  velocity  of  0.6  m  /  s  at  the  peak  of  the  winter  monsoon.  The 
simulated  thermohaline  structure  for  summer  and  winter  is  nearly  horizontal  from  east  to  west  except  at 
the  coastal  regions.  Coastal  upwelling  and  downwelling  are  also  simulated:  localized  lifting  (descending) 
of  the  isotherms  and  isohalines  during  summer  (winter)  at  the  west  boundary.  The  simulation  is  reason¬ 
able  comparing  to  the  observations.  Sensitivity  experiments  were  designed  to  investigate  the  driving 
mechanisms.  Non-linearity  is  shown  to  be  important  to  the  transport  of  baroclinic  eddy  features,  but 
otherwise  insignificant.  Transport  from  lateral  boundaries  is  of  considerable  importance  to  summer  cir¬ 
culation  and  thermal  structure,  with  lesser  effect  on  winter  monsoon  hydrology.  In  general,  seasonal  cir¬ 
culation  patterns  and  upwelling  phenomena  are  determined  and  forced  by  the  wind,  while  the  lateral 
boundary  forcing  plays  a  secondary  role  in  determining  the  magnitude  of  the  circulation  velocities. 
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1.  INTRODUCTION 

The  South  China  Sea  (SCS)  is  a  semi-enclosed  tropical  sea  located  between  the  Asian 
land  mass  to  the  north  and  west,  the  Philippine  Islands  to  the  east,  Borneo  to  the  southeast, 
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deep  basin  is  1900  km  along  its  major  axis  (northeast-southwest)  and  approximately  1 100  km 
along  its  minor  axis,  and  extends  to  over  4000  m  deep.  The  Sunda  Shelf  is  the  submerged 
connection  between  southeast  Asia,  Malaysia,  Sumatra,  Java,  and  Borneo  and  is  100  m  deep 
in  the  middle;  the  center  of  the  Gulf  of  Thailand  is  about  70  m  deep. 

The  SCS  is  subjected  to  a  seasonal  monsoon  system121.  From  April  to  August,  the  weaker 
southwesterly  summer  monsoon  winds  result  in  a  wind  stress  of  over  0.1  N/m2(Fig.  2a) 
which  drives  a  northward  coastal  jet  off  Vietnam  and  anticyclonic  circulation  in  the  SCS  (Fig. 


3a).  From  November  to  March,  the  stronger  northeasterly  winter  monsoon  winds  corre¬ 
sponds  to  a  maximum  wind  stress  of  nearly  0.3  N  /  m2  (Fig.  2b)  causing  a  southward  coastal 
jet  and  cyclonic  circulation  in  the  SCS  (Fig.  3b).  The  transitional  periods  are  marked  by  high¬ 
ly  variable  winds  and  surface  currents. 
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FIGURE  2.  Climatological  wind  stress  (0.1  N/  m2)  for  (a)  June,  and  (b)  December 
(after  Ref.  (3)). 


The  observed  circulation  patterns  of  the  intermediate  to  upper  layers  of  the  SCS  are 
primarily  forced  by  the  local  monsoon  systems121,  with  contributions  from  the  Kuroshio  Cur¬ 
rent  via  the  Bashi  Channel,  in  the  southern  half  of  the  Bashi  Channel.  The  Kuroshio  enters 
the  SCS  through  the  southern  side  of  the  channel  then  executes  a  tight,  anticyclonic  turn  and 
exits  the  SCS  near  Taiwan.  An  estimated  8-10  Sv  (1  Sv=  106  mVs)  of  the  intrusion  passes 
through  the  Bashi  Channel1'1.  This  flow  exerts  a  strong  influence  on  the  properties  of  the 
northern  SCS  waters  and  is  believed  to  contribute  to  currents  in  the  Taiwan  Strait14*. 

Eddy  behavior  in  the  SCS  has  two  distinct  features.  First,  the  number  of  cold  eddies  is  far 
greater  than  the  number  of  warm  eddies  and  second,  the  eddies  are  significantly  affected  by 
the  bottom  topography  and  are  most  likely  to  occur  near  localized  high  current  velocity*1*. 
Small  scale  eddies  with  seasonal  dependence  have  been  found  off  coastal  Vietnam 
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in  summer,  near  Natuna  Island  and  in  reef  areas.  Large  scale  eddies  have  been  found 
primarily  during  the  summer  monsoon1".  Based  on  more  complete  data  sets  such  as  U.S.  Na¬ 
vy's  Master  Oceanographic  Observational  Data  Set  (MOODS)  and  National  Centers  for 
Environmental  Prediction  (NCEP)  monthly  sea  surface  temperature  (SST)  data  set,  Chu  et 
alt5<l  identified  multi-eddy  structure  in  the  SCS. 

Upwelling  and  downwelling  occur  off  the  coast  of  central  Vietnam  and  eastern  Hainan. 
The  summer  monsoon  is  believed  to  cause  an  Ekman-type  drift  current  and  corresponding 
offshore  transport,  leading  to  upwelling1".  Others1"  have  pointed  out  that  upwelling  also  oc¬ 
curs  in  these  areas  during  winter  monsoon  conditions,  casting  doubt  on  previous  ideas  con¬ 
cerning  the  origin  of  these  features. 

In  the  north,  the  waters  are  cold  and  saline.  The  annual  variability  of  salinity  is  small, 
due  to  the  inflow  and  diffusion  of  high  salinity  water  from  the  Pacific  Ocean  through  the 
Bashi  Channel.  In  the  south  the  tropical  conditions  cause  the  waters  to  be  warmer  and 
fresher.  The  interface  of  the  southern  high  temperature,  low  salinity  w'ater  and  the  northern 
low  temperature,  high  salinity  causes  a  definite  vertical  gradient,  usually  strongest  between 
100  and  300  m  deep.  During  the  transitions  the  central  region  is  alternately  subjected  to  high 
and  low  salinity  inflow  as  the  monsoons  reverse,  resulting  in  a  region  of  higher  horizontal 
gradient  and  annual  variability.  Mixed  layer  depths  vary  from  30  to  40  m  during  the  summer 
monsoon,  and  70  to  90  m  during  the  winter  monsoon  with  variation  due  to  both  wind  and 
current’2'. 

Three-dimensional  ocean  models  have  been  used  several  times  in  the  past  to  simulate  the 
circulation  in  the  SCS.  Integrating  a  semi-implicit  12-layer  shallow  water  model  with  50  km 
resolution  under  both  mean  winter  and  summer  conditions  for  15  days,  Pohlmann’8’  simu¬ 
lated  the  reversal  of  the  upper  layer  circulation  between  the  summer  and  winter  monsoon  sea¬ 
sons.  Integrating  a  23-level  primitive  equation  model  developed  at  the  Princeton 
University”1  with  20  km  horizontal  resolution  under  monthly  mean  climatological  wind  forc¬ 
ing’"  for  a  year,  Chu  et  al.’10’  simulated  the  seasonal  variation  of  the  SCS  circulation  and  Chu 
and  Chang’1"  explained  the  formation  of  the  SCS  warm-core  eddy  in  boreal  spring.  Inte¬ 
grating  a  global  1.5-layer  reduced  gravity  model  with  0.5°  resolution,  Metzger  and 
Hurlburt’12’  successfully  simulated  the  upper  layer  circulation  and  the  mass  exchange  between 
the  SCS,  the  Sulu  Sea,  and  the  Pacific  Ocean. 

The  objective  of  this  study  is  to  simulate  the  SCS  thermohaline  structure  as  well  as  the 
circulation  and  to  investigate  physical  processes  causing  seasonal  variability.  We  used  the 
Princeton  Ocean  Model  (POM)  to  examine  the  mechanisms  causing  seasonal  variation  of  the 
SCS  circulation  and  thermal  structure.  The  control  run  is  designed  to  best  simulate  reality 
against  which  each  experiment  is  compared.  In  the  experiments,  various  external  and  internal 
factors  are  modified  and  the  resulting  circulation  patterns  and  magnitudes  are  compared  to 
the  control  run  results.  Specifically  we  estimate  the  contribution  (in  terms  of  volume 
transport,  sea  surface  elevation  and  circulation  patterns)  of  non-linear  advection,  wind  forc¬ 
ing  and  lateral  boundary  transport  to  the  ocean  features  identified  in  the  control  results. 
From  this  we  can  estimate  the  relative  importance  of  these  factors  to  SCS  oceanography. 

2.  THE  NUMERICAL  OCEAN  MODEL 

2. 1  Model  Description 

Coastal  oceans  and  semi-enclosed  seas  are  marked  by  extremely  high  spatial  and  tempo¬ 
ral  variabilities  that  challenge  the  existing  predictive  capabilities  of  numerical  simulations. 
The  POM  is  a  time  dependent,  primitive  equation  circulation  model  on  a  three-dimensional 
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grid  that  includes  realistic  topography  and  a  free  surface11*1.  Tidal  forcing  was  not  included  in 
this  application  of  the  model,  since  high  frequency  variability  of  the  circulation  is  not  consid¬ 
ered.  River  outflow  is  also  not  included.  However,  the  seasonal  variation  in  sea  surface 
height,  temperature,  salinity,  circulation  and  transport  are  well  represented  by  the  model. 
From  a  series  of  numerical  experiments,  the  qualitative  and  quantitative  effects  of 
non-linearity,  wind  forcing  and  lateral  boundary  transport  on  the  SCS  are  analyzed,  yielding 
considerable  insight  into  the  external  factors  affecting  the  region  oceanography.  The  horizon¬ 
tal  spacing  of  0.179°  by  0.175°  (approximately  20  km  resolution)  and  23  vertical  sigma 
coordinate  levels.  The  model  domain  is  from  3.06°S  to  25.07°N,  and  98.84°E  to  121.16°E, 
which  encompasses  the  SCS  and  the  Gulf  of  Thailand,  and  uses  realistic  bathymetry  data 
from  the  Naval  Oceanographic  Office  DBDB5  database  (5  minute  by  5  minute  resolution). 
Consequently,  the  model  contains  125  x  162  x  23  horizontally  fixed  grid  points.  The  horizon¬ 
tal  diffusivities  are  modeled  using  the  Smagorinsky11’1  form  with  the  coefficient  chosen  to  be 
0.2  for  this  application.  The  bottom  stress  ih  is  assumed  to  follow  a  quadratic  law 

T/>  =  D  I  y  h  |V*,  U) 

where  p„(  =  1025  kg  /  m3)  is  the  characteristic  density  of  the  sea  water,  VA  is  the  horizontal 
component  of  the  bottom  velocity,  and  CD  is  the  drag  coefficient  which  is  specified  as 
0.0025  191  in  our  model. 

2.2  Atmospheric  Forcing 

The  atmospheric  forcing  for  the  SCS  application  of  the  POM  includes  mechanical  and 
thermohaline  forcing.  The  wind  forcing  is  depicted  by 


where  (u,  v)  and  (To.r,Tot  )  are  the  two  components  of  the  water  velocity  and  wind  stress 
vectors,  respectively.  The  wind  stress  at  each  time  step  is  interpolated  from  monthly  mean 
climate  wind  stress'31,  which  was  taken  as  the  value  at  the  middle  of  the  month.  The  wind 
stress  has  a  typical  magnitude  of  0. 1-0.2  N/  m2  (Fig.  2).  Over  the  two  monsoon  seasons  the 
wind  varies  with  location  and  time,  leading  to  a  complicated  distribution  of  wind  stress. 
Surface  thermal  forcing  is  depicted  by 

KhJz  =  (pcv")*  a2C(0OBS~  0)>  <3> 

*s7£=«.2s+«2C(Sobs-S),  (4) 

where  0Obs  and  Sobs  are  the  observed  potential  temperature  and  salinity,  cp  is  the  specific 
heat,  and  Q  H  and  Q s  are  surface  net  heat  and  salinity  fluxes,  respectively.  The  relaxation 
coefficient  C  is  the  reciprocal  of  the  restoring  time  period  for  a  unit  volume  of  water. 
The  parameters  (a | ,  a 2 )  are  (0,  l)-type  switches:  ct|  =  1,  tx2  =  0,  would  specify  only  flux  forc¬ 
ing  being  applied;  ot|  =0,  a2  =  1,  would  specify  that  only  restoring  type  forcing  is  applied.  In 
this  study,  the  surface  thermal  forcing  is  determined  solely  by  restoring  forcing,  that  is  a|  =0 
and  a2  - 1  in  Eqs.  (3)-(4).  The  relaxation  coefficient  C  is  taken  to  be  0.7  m/  day,  which  is 
equivalent  to  a  relaxation  time  of  43  days  for  an  upper  layer  30  m  thick1'41.  The  net  effect  is  to 
prevent  any  deviation  from  climatology  and  ensure  that  the  SCS  acts  as  a  heat  source. 

2.3  Lateral  Boundary  Forcing 

Closed  lateral  boundaries,  i.e.,  the  modeled  ocean  bordered  by  land,  were  defined  using  a 
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free  slip  condition  for  velocity  and  a  zero  gradient  condition  for  temperature  and  salinity.  No 
advective  or  diffusive  heal,  salt  or  velocity  fluxes  occur  through  these  boundaries. 

Open  boundaries,  where  the  numerical  grid  ends  but  the  fluid  motion  is  unrestricted, 
were  treated  as  radiative  boundaries.  Volume  transport  through  the  Luzon  Strait,  Taiwan 
Strait,  and  Gasper  /  Karimata  Strait  was  defined  according  to  observations  (Table  1).  How¬ 
ever,  the  Balabac  Channel,  Mindoro  Strait,  and  Strait  of  Malacca  are  assumed  to  have  zero 
transport.  When  the  water  flows  into  the  model  domain,  temperature  and  salinity  at  the  open 
boundary  are  likewise  prescribed  from  the  climatological  data1151.  When  water  flows  out  of  the 
domain,  the  radiation  condition  was  applied, 

£(0,S)+  t/„^(0,S)=  0,  (5) 

w-here  the  subscript  n  is  the  direction  normal  to  the  boundary. 

TABLE  1.  Bi-monthly  variation  of  volume  transport  (Sv)  at  the  lateral  open 


boundaries.  The  positive  /  negative  values  mean  outflow'  /  inflow  and 
w'ere  taken  from  Ref.[2], 


Month 


Gaspar-Karimata  Straits 


Luzon  Strait 


Taiwan  Strait 


February 


August 

October 

-3.0 

1.0 

2.5 

-0.6 

December 


-3, 


-0.9 


For  the  surface  elevation,  17,  we  use  the  radiation  condition  for  the  open  boundary, 

(6) 

dt  dn  '  K01 

where  c  is  the  local  shallow  water  wave  speed,  and  F  is  the  forcing  term  including  tides.  In  this 
study,  we  use  F= 0. 

2.4  Initial  Conditions  and  Initialization 

The  model  was  integrated  with  all  three  components  of  velocity  (u,  v,  tv)  initially  set  to 
zero,  and  temperature  and  salinity  specified  by  interpolating  climatology  data1151  to  each  mod¬ 
el  grid  point.  The  model  year  consists  of  360  days  (30  days  per  month),  day  361  corresponds 
to  1  January.  It  was  found  that  90  days  were  sufficient  for  the  model  kinetic  energy  to  reach 
quasi-steady  state  under  the  imposed  conditions  (Fig.  4).  In  order  to  first  capture  the  winter 
monsoon  to  summer  monsoon  transition,  the  model  was  started  from  day  300  (30  October), 
and  run  to  day  390  (30  January  next  year)  for  spin-up.  After  day  390,  the  model  was  run  an¬ 
other  450  days  for  each  experiment. 

2.5  Mode  Splitting 

For  computational  efficiency,  the  mode  splitting  technique191  is  applied  with  a  barotropic 
time  step  of  25  seconds,  based  on  the  Courant-Friederichs-Levy  (CFL)  computational  stabil¬ 
ity  condition  and  the  external  wave  speed;  and  a  baroclinic  time  step  of  900  seconds,  based  on 
the  CFL  condition  and  the  internal  wave  speed. 


2.6  Experiment  Design 


Our  approach  was  to  carry  out  four  numerical  experiments:  one  control  and  three  sensi¬ 
tivity  runs.  All  runs  were  completed  for  the  same  18  month  period  encompassing  both  sum¬ 
mer  and  winter  monsoons,  and,  except  as  specified  below,  utilized  the  same  initial  conditions 
(on  day  300).  Run  1  was  the  control  run.  The  three  sensitivity  runs  were:  Run  2,  linear 


***** 


3.  THE  SIMULATION  (CONTROL  RUN) 


3.1  Circulation 


FIGURE  4.  Temporal  variation  or  total  kinetic  energy.  The  model  is  integrated  from 
November  and  reaches  quasi-steady  state  after  90  days  of  integration. 
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During  the  summer  monsoon  period  (mid-May  to  August)  winds  blow  from  the 
southwest  and  the  SCS  surface  (we  refer  to  z  =  0)  circulation  generally  follows  suit  with 
anticyclonicity  in  the  southern  basin  (Figs.  5a-5c).  Inflow  is  through  the  southern  Gaspar 
and  Karimata  Straits  and  outflow  is  through  the  northern  Taiwan  Strait  and  eastern  Bashi 
Chanel  (Table  1).  The  simulated  summer  (June-August)  mean  general  circulation  pattern  has 
the  following  features.  Velocities  reach  1  m  /  s  at  the  peak  of  the  summer  monsoon  within  the 
Vietnam  Coastal  Jet  (VCJ).  The  western  boundary  current  splits  into  two  currents  at  12°N: 
the  coastal  current  and  off-shore  current.  The  off-shore  current  further  bifurcates  into 
northward  and  southeastward  branches.  The  bifurcation  point  is  along  110°E  at  14°N  in 
May,  and  southward  shifts  to  10°N  in  July.  The  cross  basin  zonal  current  reaches  14°N  with  a 
core  speed  of  0.3  m/  s  in  May  and  shifts  to  12°N  with  a  core  speed  of  0.6  m/  s  in  July.  The 
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FIGURE  5.  Surface  circulation  for  control  run  on  (a)  May  15;  (b)  July  15;  (c)  Sep¬ 
tember  15;  (d)  October  15;  (e)  December  15;  and  (c)  March  15. 
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coastal  branch  continues  north  then  east  at  Hainan  Island.  Such  a  VCJ  summer  bifurcation 
was  simulated  using  a  barotropic  model1'4171. 

During  the  winter  monsoon  period  (November  to  March)  the  winter  Asian  high  pressure 
system  brings  strong  winds  from  the  northeast  and  the  SCS  surface  circulation  pattern  is 
cyclonic  (Figs.  5d-5f).  Inflow  from  the  Bashi  Channel  (the  Kuroshio  intrusion)  and  Strait  of 
Taiwan  augments  currents  southwest  along  the  Asian  continental  shelf,  then  south  along  the 
coast  of  Vietnam  and  eventually  out  through  the  Gasper  and  Karimata  Straits  in  the  south. 
Western  intensification  of  the  general  cyclonic  circulation  pattern  was  also  simulated.  From 
the  south  coast  of  Hainan  Island,  the  current  intensifies  as  it  flows  from  north  to  south  along 
the  Vietnam  coast.  Average  speed  is  around  0.8  m/s  in  the  core.  As  winter  progresses,  near 
surface  currents  along  the  coast  of  Borneo  begin  to  turn  northwest,  eventually  flowing 
directly  away  from  the  coastline  and  into  the  Natuna  eddy  where  the  southern  edge  of  the 


115  120 

Longitude  (E)  Longitude  (E) 

FIGURE  6.  Mean  surface  circulation  for  control  run  during  (a)  summer,  and  (b) 
winter. 


deep  basin  meets  the  Sunda  Shelf  (Fig.  5e);  the  south  extension  of  the  VCJ  veers  west.  In  early 
spring  a  northward  surface  current  developed  adjacent  to  the  VCJ  (Fig.  5f)  and  persists 
through  the  transition  periods  to  the  start  of  the  summer  monsoon. 

In  order  to  compare  the  difference  of  the  two  monsoon  seasons,  we  averaged  the  model 
variables  such  as  velocity,  surface  elevation,  temperature,  and  salinity  over  the  two  periods: 
June-August  (summer)  and  December-February  (winter.)  The  surface  circulations  for  the 
summer  (Fig.  6a)  show  an  anticyclonic  gyre  in  the  deep  basin  with  the  VCJ  flowing  northward 
to  northeastward.  The  winter  circulations  (Fig.  6b)  indicate  a  reverse  pattern:  a  cyclonic  gyre 
in  the  deep  basin  with  the  VCJ  flowing  southward  to  southeastward.  The  seasonal  variation 
of  VCJ  is  illustrated  by  the  zonal  cross-sections  of  v-component  at  13°N  for  winter  (Decem¬ 
ber  -February)  and  summer  (June  -  August).  The  summer  VCJ  is  flowing  northward  with  a 
mean  maximum  speed  of  0.5  m  /  s  and  a  width  of  100  km,  extending  to  a  depth  of  around  200 
m  (Fig.  7a).  The  winter  VCJ  is  a  much  stronger  southward-flowing  boundary  current  with  a 
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FIGURE  7.  Mean  latitudinal  velocity  for  control  run  at  the  13°N  cross-section  dur¬ 
ing  (a)  winter,  and  (b)  summer. 
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FIGURE  8.  Mean  surface  elevation  (m)  for  control  run  during  (a)  summer,  and  (b) 
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mean  maximum  velocity  of  0.95  m  /  s  and  a  width  of  100  km.  extending  to  a  depth  of  500  m 
(Fig.  7b). 

Sea  surface  heights  for  the  winter  show  a  0.2  m  depression  over  the  deep  basin  with  posi¬ 
tive  height  in  the  south  over  the  Sunda  Shelf  and  the  Gulf  of  Thailand  (Fig.  8a).  Average 
summer  sea  surface  elevation  varies  from  -0.1  to  0.1  m  with  a  southeast  upward  till  (Fig.  8b). 


3.2  Thermohaline  Structure 


Isotherms  and  isohalines  for  summer  and  winter  are  nearly  horizontal  from  east  to  west 
except  at  the  coastal  regions  (Fig.  9).  Localized  lifting  of  the  isotherms  and  isohalines  during 
summer  at  the  west  boundary  indicates  coastal  upwelling  inside  the  VCJ  (Figs.  9a  and  9c). 
Localized  descending  of  the  isotherms  and  isohalines  during  winter  at  the  west  boundary  in¬ 
dicates  coastal  downwelling  inside  the  VCJ  (Figs.  9b  and  9d).  Both  coastal  upwelling  and 
downwelling,  combined  with  the  high  velocity  shear  across  the  VCJ,  result  in  baroclinic  insta¬ 
bility. 
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FIGURE  9.  Thermohaline  feature  at  the  I3°N  zonal  cross-section  simulated  by  the 
control  run:  (a)  summer  temperature  (°C);  (b)  winter  temperature  (°C);  (c)  summer 
salinity  (psu),  and  (d)  winter  salinity  (psu). 


Latitudinal  thermohaline  variation  for  summer  and  winter  is  shown  in  Fig,  10.  In  general 
model  thermohaline  structure  is  consistent  with  the  two  SCS  water  masses  described  by 
Wyrtki121.  Over  the  southern  basin  there  is  a  general  lifting  of  isotherms  of  40-50  m  from  win¬ 
ter  to  summer  above  200  m.  In  the  northern  SCS,  near-surface  waters  (20-50  m)  are  influ- 
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cnced  by  the  winter  inflow'  of  North  Pacific  Kuroshio  water.  They  are  uniformly  colder  and 
more  saline  and  there  is  a  weaker  thermocline  in  winter  than  in  summer.  In  the  south  the 
equatorial  climate  and  summer  inflow  from  the  shallow  Java  Sea  cause  the  water  mass  to  be 
fresher,  warmer  with  a  slightly  deeper  thermocline.  The  model  SCS  summer  thermocline 
depth  is  in  good  agreement  with  Ref.  [2],  at  30  and  50  m,  as  opposed  to  the  model  winter 
depth  which  is  much  shallower. 
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FIGURE  10.  Same  as  Fig. 9  except  for  the  1 13°E  latitudinal  cross-section. 
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3.3  Volume  Transport 


The  volume  transport  stream  functions  for  summer  and  winter  are  shown  in  Figs.  11a 
and  1  lb.  For  the  SCS  deep  basin,  the  barotropic  flow  reveals  an  anticyclonic  (cyclonic)  gyre 
in  the  summer  (winter).  For  the  SCS  continental  shelf/ slope  regions,  the  volume  transport 
has  an  evident  seasonal  variation.  Winter  volume  transport  is  strongest  along  the  north  and 
west  slope  with  significant  eddy  and  meander  activity,  and  across  the  southern  slope  of  the 
deep  basin.  Average  winter  VCJ  transport  is  about  10.6  Sv.  Summer  shows  a  weaker  VCJ 
transport  (5.5  Sv)  with  similar  transport  along  the  southern  slope  of  the  deep  basin,  again 
with  noticeable  eddy  activity  (Table  2). 


3.4  Natuna  Island  Eddy 


During  winter,  the  confluence  of  the  VCJ  and  the  offshore  flow  from  Borneo  in  the 
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4.  DRIVING  MECHANISMS 


We  analyzed  the  results  of  the  four  experiments  to  identify  the  driving  mechanisms  for 
the  SCS  circulation. 


494 


C  hinese  Journal  of  Atmospheric  Sciences 


TABLE  2.  Volume  transport  (Sv)  of  individual  features 


simulated  by  the  control  run. 


_ Feature 

Summer  Vietnam  Coastal  Jet 
Winter  Vietnam  Coastal  Jet 
Natuna  Island  Eddy _ 


Transport 
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106 
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FIGURE  11.  Mean  total  volume  transport  streamfunction  (Sv)  for  control  run  dur¬ 
ing  (a)  summer,  and  (b)  winter. 


Sunda  Shelf  generates  a  mesoscale  cyclonic  eddy  (Fig.  6b),  here  termed  the  Natuna  Island 
Eddy  (NIE),  consistent  with  observations  of  late  winter  surface  circulation  described  by 
Wyrlki  (Fig.  3b).  This  eddy  shows  little  variability  during  the  winter,  with  maximum  swirl  ve¬ 
locity  of  0.6  m/  s  at  the  peak  of  the  winter  monsoon  but  no  variation  in  position.  NIE  is  a 
weak  anticyclonic  eddy  (Fig.  6a)  during  summer,  and  a  strong  cyclonic  eddy  (Fig.  6b)  during 
winter.  The  summer  NIE  has  an  average  core  velocity  of  0. 1  m  /  s.  The  winter  NIE  has  an  av- 
erage  core  velocity  of  0.45— 0.5  m/s.  The  east— west  slanted  eddy  core  is  deeper  on  the  west 
side  (depth  ~  40  m)  of  the  eddy  than  on  the  east  side  (depth  ~  20  m).  The  vertical  shear  of 
horizontal  velocity  promotes  baroclinic  instability. 


4. 1  Effects  of  Non-Linearity  (  Run  1  -  Run  2) 

In  the  first  sensitivity  study  the  nonlinear  advection  terms  were  removed  from  the 
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FIGURE  12.  Mean  latitudinal  velocity  for  control  run  at  the  6°N  cross-section  dur¬ 
ing  (a)  summer,  and  (b)  winter. 
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FIGURE  13.  Mean  surface  circulation  for  linear  run  during  (a)  summer,  and  (b)  win 
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dynamic  equations.  Otherwise  the  same  parameters  as  the  control  run  were  used.  The  near 
surface  vector  velocities  for  the  summer  (Fig.  13a)  and  winter  (Fig.  13b)  illustrate  the  similari¬ 
ty  to  the  control  run  results  (Fig.  6),  which  indicates  that  the  non-linearity  does  not  change 
the  general  circulation  pattern.  However,  our  computation  shows  that  the  non-linearity 
causes  a  noticeable  change  in  the  volume  transport  and  the  strength  of  the  western  boundary 


currents. 

Figs.  14a  and  15a  are  plots  of  the  difference  in  the  volume  transport  streamfunction  ( A«/r ) 
for  the  control  run  minus  the  linear  run  and  represent  nonlinear  effects  for  summer  and  win¬ 
ter,  respectively.  The  biggest  Ai/>  during  both  seasons  is  a  deep-basin  double-gyre  feature:  an 
anticyclonic  gyre  (-5  Sv)  west  of  Luzon  and  a  cyclonic  gyre  (10  Sv)  west  of  Liyue  Bank  (1 1°N, 

1 15°E).  Such  a  double-gyre  structure  indicates  (a)  a  westward  cross-basin  transport  (15  Sv) 
near  13°N  (between  the  two  gyres),  (b)  increase  (decrease)  of  cyclonic  (anticyclonic)  transport 
in  the  southwestern  basin  (south  of  13°N),  and  (c)  increase  (decrease)  of  anticyclonic 
(cyclonic)  transport  in  the  northeastern  basin  (north  of  13°N).  Furthermore,  two  nearly  paral¬ 
lel  zero  curves  of  Aip  show  up  near  the  Vietnam  coast  during  the  summer  (Fig.  14a).  One 
curve  is  close  to  the  coast,  and  the  other  curve  is  off  shore.  This  shows  that  Aip  =0  on  the  en¬ 
tire  Vietnam  shelf,  resulting  in  no  change  in  the  northward  VCJ  transport  in  the  summer. 
While,  during  the  winter  such  a  feature  (At/'  =0)  appears  only  north  of  13°N  and  some  south¬ 
ward  volume  transport  increment  (At//  =4  Sv)  shows  up  near  the  Vietnam  coast  south  of 
13°N.  After  averaging  At/t  over  the  Vietnam  coast  shelf  break,  we  estimated  that  the  nonlinear 
advection  brings  a  3.1  Svof  the  VCJ. 

The  evident  nonlinear  effect  on  the  western  boundary  current  can  be  identified  by  the  ve¬ 
locity  difference  for  the  control  run  minus  the  linear  run.  A  cross-section  of  the  latitudinal  ve¬ 
locity  difference  during  the  summer  at  6°N  shows  that  nonlinear  advection  augments  coastal 
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FIGURE  14.  Summer  volume  transport  anomaly  (Sv)  caused  by  effects  of  (a) 
nonlinearity;  (b)  winds;  and  (c)  boundary  forcing. 
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core  current  at  the  western  boundary  by  about  0.25  m/  s  and  is  responsible  for  offshore  eddy 
activity  (Fig.  16a).  Such  a  nonlinear  effect  was  also  found  in  the  surface  velocity  vector  differ¬ 
ence  (AV)  field.  The  primary  feature  of  the  nonlinear  effect  on  circulation  during  the  winter 
(Fig.  17a)  is  similar  to  that  during  the  summer  (Fig.  18a):  quite  evident  near  the  west  coast 
and  not  evident  elsewhere. 

4.2  Wind-Induced  Circulation  (  Run  I  -Run  3 ) 

The  difference  between  Run  1  and  Run  3  shows  the  monsoon  wind  effects  on  the  SCS 
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FIGURE  16.  Summer  latitudinal  velocity  anomaly  (m/s)  at  the  6°N  cross-section, 
caused  by  effects  of  (a)  nonlinearity,  (b)  winds,  and  (c)  boundary  forcing. 


105  110  115  120  100 

Longitude  (E) 

105  no  115  12< 

Longitude  (E) 

> 

100  105  110  115 

Longitude  (E) 

FIGURE  15.  Winter  volume  transport  anomaly  (Sv)  caused  by  effects  of  (a) 
nonlinearity;  (b)  winds:  and  (c)  boundary  forcing. 
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FIGURE  17.  Summer  surface  current  anomaly  caused  by  efTects  of  (a)  nonlinearity, 
(b)  winds,  and  (c)  boundary  forcing. 
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FIGURE  18.  Winter  surface  current  anomaly  caused  by  efTects  of  (a)  nonlinearity 
(b)  winds,  and  (c)  boundary  forcing. 

surface  circulation  and  SST  for  the  summer  and  winter  respectively.  The  primary  ft 
the  summer  wind  effect  (Fig.  17b)  on  circulation  and  volume  transport  is  in  contra 


499 


flows  with  0.3— 0.4  m/s  velocity  (Fig.  18b)  corresponding  to  a  southward  volume  transport 
of  4.6  Sv  (Fig.  15b).  The  China  slope  flow  is  actually  slightly  more  than  the  control  run  but 
this  may  be  accounted  for  by  nonlinearity.  When  winds  are  removed  from  the  forcing,  the 
model  results  show  considerable  change  in  the  structure  of  the  summer  circulation  and  a  vir¬ 
tual  disappearance  of  features  in  the  winter  circulation  (note  similarities  in  Figs.  6b  and  18b); 
surface  velocities  are  reduced  by  0.20-0.3  m/  s. 

The  summer  w'inds  drive  a  mesoscale  anticyclonic  eddy  near  a  shallow  depression  in  the 
Sunda  shelf  and  Natuna  Island,  with  rotation  speeds  of  0.15-0.2  m/  s  (Fig.  16b.)  Winds  also 
induce  the  loop  current  joining  the  Sunda  shelf  flow  and  the  VCJ,  bifurcation  of  the  coastal 
current,  and  cross-basin  flow  (Fig  17b.)  The  corresponding  wind  effect  on  summer  volume 
transport  is  largely  limited  to  the  northern  SCS  (Fig.  14b)  with  a  6-8  Sv  maximum  transport 
along  the  China  slope  and  a  4-6  Sv  eastward  transport  centered  on  13°N,  a  region  of  signifi¬ 
cant  nonlinear  effect. 

Summer  wind  reduces  SST  slightly  (1-2°C  )  from  upwelling  off  central  and  northern 
Vietnam  and  the  east  coast  of  Hainan  (Fig.  19a).  Winter  wind  decreases  surface  temperature 
to  a  lesser  extent  on  the  east  side  of  the  basin  along  the  Palawan  Trough  (Fig.  19b). 
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FIGURE  19.  Sea  surface  temperature  anomaly  caused  by  wind  effect. 


120 


Volume  22,  No.4 


Summer  winds  shift  the  average  surface  elevation  to  the  north,  increasing  the  elevation 
5-10  cm,  while  elsewhere  decreasing  the  height  by  the  same  magnitude  (Fig.  20b).  The  winter 
wind  increases  the  surface  elevation  over  the  continental  shelves  by  0.15-0.18  m,  and  de¬ 
creases  the  elevation  over  the  deep  basin  (Fig.  2lb). 


4.3  Inflow/  Outflow  Induced  Circulation  (  Run  I  -  Run  4) 

The  third  sensitivity  study  used  the  control  run  equations  and  forcing  but  closed  all  open 
lateral  boundaries,  preventing  transport  of  mass,  heat  or  salinity  through  the  Luzon  Strait, 
Taiwan  Strait  or  Gasper  and  Karimata  Straits.  With  no  inflow  or  outflow  the  summer 
anticyclonic  gyre  and  winter  cyclonic  gyre  are  more  pronounced.  Increased  recirculation  gen¬ 
erally  leads  to  greater  horizontal  and  vertical  variability  of  the  current  structure. 


(a) 

Summer  SST  (Cntrt-No  Wind) 


(b) 

Winter  SST  (Cntrt-No  Wind) 
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The  results  for  summer  lateral  boundary-forcing  retain  more  of  the  anticyclonic  nature 
of  the  circulation  than  do  those  from  wind-forcing,  see  volume  transport  streamfunction  dif¬ 
ference  Ai/f  (Figs.  14b  and  14c.)  The  flow  is  generally  slope  hugging.  The  largest  gradient  of 
Ai p  is  located  off  South  Vietnam  (Fig.  14c),  causing  the  strong  western  boundary  currents. 
The  southern  inflow  splits  south  of  Natuna  Island  and  fans  across  the  Sunda  shelf,  forming  a 
small  anticylonic  eddy  collocated  with  that  forced  by  the  summer  wind.  The  VCJ  transports 
5.5  Sv  and  is  accompanied  by  a  counter  current  nearby.  At  6°N  (Fig.  16c),  however,  notice 
that  the  average  northward  velocity  throughout  most  of  the  shallow  shelf  region  shows  no 
discernible  pattern  of  change,  while  the  shelf  current  along  the  100  m  isobath  decreases  by 
0.05-0.1  m/  s.  The  summer  circulation  patterns  change  dramatically  when  the  lateral  bound¬ 
ary'  transport  is  removed.  This  can  be  seen  from  the  surface  velocity  vector  difference 
AV  (Fig.  17c).  There  is  still  western  intensification  of  current  along  the  coast  of  Malaysia,  and 
it  still  joins  the  flow  out  of  the  Gulf  of  Thailand  to  contribute  to  an  intensification  of  current 
off  southern  Vietnam. 

(a)  0>)  (0 


Summer  Elev  (CnM-Unear)  Summer  Elev  (Cntrt-No  Wind)  Summer  Elev  (Cmn-Clsd  Bdys) 


100  105  110  IIS  120  100  105  110  IIS  120  100  105  110  IIS  120 


Longitude  [E]  Longitude  [E]  Longitude  [E] 

FIGURE  20.  Summer  surface  elevation  anomaly  (m)  caused  by  effects  of  (a)  winds, 
and  (b)  boundary  forcing. 

Winter  closed  boundary  circulation  patterns  on  the  other  hand,  show  less  difference  in 
structure  from  the  control  run  but  more  variability  in  magnitude.  There  is  an  average  decrease 
of  0.1-0.4  m/s  in  current  speed  of  the  winter  VCJ  (Fig.  15c),  which  equates  to  4-6  Sv 
southward  volume  transport.  The  positive  wind  curl  acts  to  strengthen  the  cyclonic  gyre  na¬ 
ture  of  overall  circulation.  The  Kuroshio  intrusion  and  inflow  through  the  Taiwan  Strait  are 
obviously  supplemented  by  recirculation  flow  from  along  the  coast  of  Luzon  Island.  Near 
Natuna  Island  it  is  noteworthy  that  the  structure  of  cross-basin  circulation  and  current  flow 
away  from  the  Borneo  coast  are  unchanged.  Similarly,  it  is  apparent  that  the  spatial  extent 
and  shape  of  the  gyre  northeast  of  Natuna  Island  is  unchanged  by  closing  the  boundary  flow. 
In  cross  section,  however,  it  can  be  seen  that  the  NIE  does  lose  some  of  its  velocity  — there  is 
a  0. 1-0.3  m  /  s  decrease  in  the  average  core  velocity  on  the  western  side  at  6°N  in  this  run  cor¬ 
responding  to  a  volume  transport  of  3  Sv,  associated  with  the  decrease  in  velocity  of  the  VCJ. 
The  difference  in  velocity  is  much  smaller  on  the  eastern  side,  suggesting  that  some  other  ef- 
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feet  dominates  the  current  structure  in  that  region. 

The  nonlinearity  increases  the  surface  elevation  slightly  (0.05  m)  near  the  south  Vietnam 
coast  in  summer  (Fig.  20a),  and  has  almost  no  effect  on  the  surface  elevation  in  winter  (Fig. 
21a).  The  summer  lateral  boundary  transport  decreases  the  average  surface  elevation  over  the 
China-Vietnam  continental  shelf  by  0.05-0.15  m,  while  elsewhere  increases  the  height  with  a 
maximum  accretion  of  0.2-0.4  m  near  the  Karimata  Strait  (Fig.  20c).  The  winter  lateral 
boundary  transport  increases  the  average  surface  elevation  over  the  China-Vietnam  conti¬ 
nental  shelf  by  0.05-0.2  m,  while  elsewhere  decreases  the  height  with  a  maximum  reduction  of 
0.25  m  near  the  Karimata  Strait  (Fig.  21c). 


(6) 


(c) 
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FIGURE  21.  Winter  surface  elevation  anomaly  (m)  caused  by  effects  of  (a)  winds, 
and  (b)  boundary  forcing. 


5.  CONCLUSIONS 


(1)  The  SCS  circulation  and  thermohaline  structure  were  simulated  in  this  study  by  the 
POM  model  under  the  climatological  forcing.  During  summer  (winter)  monsoon  period  the 
SCS  surface  circulation  is  generally  anticyclonic  (cyclonic)  with  a  strong  western  boundary 
current  (width  around  100  km)  -  Vietnam  Coastal  Jet.  This  jet  has  a  strong  seasonal  variation 
and  flows  northward  (southward)  during  summer  (winter)  with  a  mean  maximum  speed  of  0.5 
m/  s  (0.95  m/  s),  mean  volume  transport  of  5.5  Sv  (10.6  Sv),  extending  to  a  depth  of  around 
200  m  (500  m).  During  summer,  the  western  boundary  current  splits  and  partially  leaves  the 
coast;  the  bifurcation  point  is  at  14°N  in  May,  and  shifts  south  to  10°N  in  July.  The  cross-ba¬ 
sin  zonal  current  reaches  12°N  and  a  core  speed  of  0.6  m/  s  by  the  end  of  the  summer.  The 
coastal  branch  continues  north  then  east  at  Hainan  Island.  Average  summer  sea  surface  eleva¬ 
tion  varies  from  -0.1  to  0.1  m  with  a  southeast  upward  tilt.  Sea  surface  heights  for  the  winter 
show  a  0.2  m  depression  over  the  deep  basin  with  positive  height  in  the  south  over  the  Sunda 
Shelf  and  the  Gulf  of  Thailand.  Besides,  the  POM  model  successfully  simulates  a  mesoscale 
eddy  in  the  Sunda  Shelf,  that  is  the  Natuna  Island  Eddy.  This  eddy  is  cyclonic  (anticyclonic) 
with  maximum  swirl  velocity  of  0.6  m  /  s  at  the  peak  of  the  winter  monsoon. 
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Isotherms  and  isohalines  for  summer  and  winter  are  nearly  horizontal  from  east  to  west 
except  at  the  coastal  regions.  Coastal  upwetling  and  downwelling  are  also  simulated:  localized 
lifting  (descending)  of  the  isotherms  and  isohalines  during  summer  (winter)  at  the  west  bound¬ 
ary.  Both  coastal  upwelling  and  downwelling  (causing  horizontal  thermohaline  gradient), 
combined  with  the  high  velocity  shear  across  the  coastal  jet,  result  in  baroclinic  instability. 
This  mechanism  may  contribute  to  the  summer  jet  bifurcation.  In  general  model  thermohaline 
structure  is  consistent  with  the  two  SCS  water  masses  described  by  Wyrtki121.  Over  the  south¬ 
ern  basin  there  is  a  general  lifting  of  isotherms  of  40-50  m  from  winter  to  summer  above  200 
m.  In  the  northern  SCS,  near-surface  waters  (20-50  m)  are  influenced  by  the  winter  inflow  of 
North  Pacific  Kuroshio  water. 

(2)  The  model  nonlinear  effects  on  circulation  in  the  summer  SCS  are  localized.  There  is 
a  dipole  in  volume  transport  over  the  western  corner  of  the  southern  slope  of  the  deep  basin  in 
conjunction  with  anticyclonic  transport  extending  across  the  mid-basin  above  this  southern 
slope.  The  winter  mid-basin  transport  is  similarly  anticyclonic  but  with  the  dipole  reducing  to 
a  single  cyclonic  center  and  the  addition  of  another  cyclonic  transport  center  near  Natuna  Is¬ 
land.  The  mid-basin  transport  and  winter  cyclonic  center  have  a  weak  surface  current  signa¬ 
ture.  The  winter  surface  current  for  nonlinear  dynamic  effects  clearly  shows  the  NIE.  Eleva¬ 
tion  is  little  effected  by  nonlinear  dynamics. 

(3)  The  model  wind  effects  on  the  SCS  are  more  perversive  than  those  of  non-linear  dy¬ 
namic  effects.  Similar  to  nonlinear  effects  there  are  areas  of  anticylonic  volume  transport 
mid-basin  in  the  summer  and  winter,  although  the  mid-basin  surface  currents  reverse.  Sum¬ 
mer  wind  effects  on  transport  are  generally  divided  along  the  axis  of  the  monsoon  wind,  SW 
to  NE,  with  anticyclonicity  and  positive  elevation  in  the  southeast  and  cyclonic  transport  and 
negative  elevation  in  the  northwest.  Winter  wind  effect  is  largely  cyclonic,  including  a  center 
(with  a  strong  surface  current  signature)  for  the  NIE,  with  a  matching  steeper  depression  in 
surface  elevation. 

(4)  The  model  boundary-forcing  effect  on  volume  transport  is  most  easily  distinguished 
as  anticyclonic  in  the  summer  and  cyclonic  in  the  winter.  The  stronger  transport  centers  of 
both  seasons  have  strong  surface  current  signatures  including  the  NIE.  Topographically  li¬ 
nked  eddy  flow  and  transport  are  evident  in  the  summer  while  winter  flow  narrowly  follows 
the  200  m  isobath  until  it  encounters  the  Sunda  Shelf  and  Natuna  Island. 

In  the  summer  SCS  volume  transport  is  mostly  due  to  boundary-forcing.  Surface 
movement  is  marginally  dominated  by  wind-forcing,  when  the  boundary-forcing  surface  ef¬ 
fect  (away  from  the  southern  boundary)  coalesces  with  associated  volume  transport.  Winds 
drive  the  offshore  summer  bifurcation  of  the  VCJ.  The  anticylonic  gyre  south  of  Hainan  is  a 
surface  phenomenon  driven  by  boundary-flow  interaction  with  topography  and  nonlinear  dy¬ 
namics. 

In  the  winter  surface  movement  due  to  boundary-forcing  is  largely  limited  to  the  perime¬ 
ter  of  the  SCS  basin,  specifically  a  narrow  swath  over  the  200  m  isobath,  away  from  the  south¬ 
ern  boundary.  Surface  flow  and  volume  transport  are  dominated  by  wind  effect  except  at  the 
NIE,  a  mesoscale  feature  that  must  be  qualified  by  nonlinear  dynamics.  Also  to  be  qualified 
by  nonlinear  dynamics  is  cross-basin  transport  for  both  seasons.  Further  consideration  of 
cross-basin  transport  should  accommodate  flow  through  the  Mindoro  Strait.  Inspection  of 
observational  surface  currents  from  the  Sulu  Sea  shows  that  these  flows  may  be  important  to 
formation  of  the  two  gyre  system  that  exists  in  the  spring. 

(5)  Future  studies  should  concentrate  on  less  simplistic  scenarios.  Realistic  surface  heat 
and  salt  fluxes  should  be  included  and  the  use  of  extrapolated  climatological  winds  needs  to 
be  upgraded  to  incorporate  synoptic  winds  to  improve  realism.  Finally,  the  assumption  of 
quasi-linearity  that  allowed  us  to  use  simple  difference  to  quantify  the  effect  of  external  fore- 
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ing  needs  to  be  rigorously  tested.  It  is  very  important  to  develop  a  thorough  methodology  to 
perform  sensitivity  studies  under  the  highly  non-linear  conditions  that  may  exist  in  the  littor¬ 
al  environment. 
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